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: ABSTRACT 
> ■ 

' The outermost layers of some neutron stars are likely to be dominated by hydrogen, 

. as a result of fast gravitational settling of heavier elements. These layers directly mediate 

thermal radiation from the stars, and determine the characteristics of X-ray/EUV spectra. 
For a neutron star with surface temperature T <^ 10^ K and magnetic field B ^ 10^^ G, 
i various forms of hydrogen can be present in the envelope, including atom, poly-molecules, 

^ I and condensed metal. We study the physical properties of different hydrogen phases on 

• the surface of a strongly magnetized neutron star for a wide range of field strength B and 

surface temperature T. Depending on the values of B and T, the outer envelope can be 
O ' either in a nondegenerate gaseous phase or in a degenerate metallic phase. For T ^ 10^ K 

^ . and moderately strong magnetic field, B ^ 10^^ G, the envelope is nondegenerate and the 

^ I surface material gradually transforms into a degenerate Coulomb plasma as density increases. 

^ • For higher field strength, B » 10^^ G, there exists a first-order phase transition from 

• i-H . the nondegenerate gaseous phase to the condensed metallic phase. The column density of 

^ ' saturated vapor above the metallic hydrogen decreases rapidly as the magnetic field increases 

■ or /and temperature decreases. Thus the thermal radiation can directly emerge from the 

degenerate metallic hydrogen surface. The characteristics of surface X-ray/EUV emission 
for different phases are discussed. A separate study concerning the possibility of magnetic 
field induced nuclear fusion of hydrogen on the neutron star surface is also presented. 

Subject headings: stars: neutron - stars: atmospheres - magnetic fields - atomic pro- 
cesses - equation of state - radiation mechanisms: thermal 
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1. INTRODUCTION 



It has long been realized that neutron stars should remain detectable as soft X-ray 
sources for ~ 10^ years after their birth (Chiu &; Salpeter 1964; Tsuruta 1964). The last 
thirty years have seen significant progress in our understanding of various physical processes 
responsible to neutron star cooling (see e.g., Pethick 1992; Umeda, Tsuruta & Nomoto 1994; 
Reisenegger 1995). The advent of imaging X-ray telescopes has now made it possible to 
observe isolated neutron stars directly by their surface radiation. In particular, ROSAT has 
detected pulsed X-ray thermal emission from a number of radio pulsars (see Becker 1995 for 
a review). Several nearby pulsars have also been detected by EUVE (Edelstein &; Bowyer 
1996). On the other hand, old isolated neutron stars (10^ — 10^ of which are thought to exist 
in the Galaxy), heated through accretion from interstellar material, are also expected to be 
common sources of soft X-ray/EUV emission (e.g, Treves & Colpi 1991; Blaes & Madau 
1993. See Walter et al. 1996 for a possible detection and Pavlov et al. 1996 for spectral 
interpretations) . It has been suggested that some of the unidentified sources detected in the 
EUVE and ROASAT WFC all-sky surveys may be associated with such old neutron stars 
(Shemi 1995). Overall, the detections of the surface emission from neutron stars have the 
potential of constraining the nuclear equation of state, various heating/accretion processes, 
magnetic field structure and surface chemical composition. Future observations are likely to 
extend to surface temperatures as low as 10^ K. Confrontation with theory requires detailed 
understanding of the physical properties of the outer layers of neutron stars, in the presence 
of intense magnetic fields {B ^ 10^^ G) and low temperatures. 

For young radio pulsars that have not accreted much gas, one might expect the surfaces 
to consist mainly of iron-peak elements formed at the neutron star birth. However, once 
the neutron stars accrete material, or have gone through a phase of accretion, either from 
the interstellar medium or from a binary companion, a hydrogen envelope will form on the 
top of the surface unless it is completely burnt out. It is this envelope that mediates the 
radiation from the neutron star surface. While the strong magnetic field and/or rapid stellar 
spin may prevent large-scale accretion, it should be noted that even with a low accretion rate 
of 10^° g s~^ (typical of accretion from the interstellar medium) for one year, the accreted 
material will be more than enough to completely shield the iron surface of a neutron star. 
The lightest elements, H and He, are likely to be the most important chemical species in the 
envelope due to their predominance in the accreting gas and also due to quick separation of 
light and heavy elements in the gravitational field of the neutron star (e.g., the settling time 
of C in a 10^ K hydrogen photosphere is of order a second). If the present accretion rate is 
low, gravitational settling produces a pure H atmosphere. If the temperature on the neutron 
star surface is not too high, light atoms, molecules and metal grains may form. The purpose 
of this paper is to study the phase diagram and the physical properties of the hydrogen 
envelope for a wide range of magnetic field strength and surface temperature. 

A strong magnetic field can dramatically change the structure of atoms, molecules and 
condensed matter (see Ruderman 1974 for an early review, and Ruder et al. 1994 for a 
recent text. Heavier atoms arc also discussed by Lieb et al. 1994). The atomic unit Bo for 
the magnetic field strength and a dimensionless parameter b are 

b=^; Bo = = 2.351 X 10^ G. (1.1) 

Bo h 

When b » 1, the cyclotron energy of the electron hUf, = U{eB / mf,c) = 11.58Si2 keV, where 
Bi2 = b/425.4 is the magnetic field strength in units of 10^^ G, is much larger than the typical 
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Coulomb energy, thus the Coulomb forces act as a perturbation to the magnetic forces on 
the electrons, and at most temperatures the electrons settle into the ground Landau level. 
Because of the extreme confinement of electrons in the transverse direction, the Coulomb 
force becomes much more effective for binding electrons in the parallel direction. The atom 
has a cigar-like structure. Moreover, it is possible for these elongated atoms to form molecular 
chains by covalent bonding along the field direction. In two recent papers (Lai, Salpeter & 
Shapiro 1992; Lai & Salpeter 1996; hereafter referred as Paper I and Paper II), we have 
studied the electronic structure and energy levels of various forms of hydrogen in strong 
magnetic field, including atoms, poly-molecular chains Hat and condensed metal. In contrast 
to Fe chains, which is unbound at zero-pressure (Jones 1985; Neuhauser, Koonin & Langanke 
1987), we found that for typical magnetic field strength, B12 ^ 1, the infinite H chains (thus 
the metallic hydrogen) are bound relative to individual atoms, and the cohesive energy 
increases with the field strength. This gives rise to the possibility of condensation of metallic 
hydrogen for sufficiently low temperature and/or high magnetic field. We quantify this phase 
transition in §4 of this paper. 

In this study, we shall focus on the magnetic field strength in the range of 10^^ ^ 

B ^ 10^^ G so that b » 1 is well-satisfied. While field strengths of order 10^^ _ iqis q 
are considered typical for most neutron stars (with the exception of old millisecond pulsars 
and neutron stars in low-mass X-ray binaries), it should be noted that the only physical 
upper limit to the neutron star magnetic field strength is the virial equilibrium value, ~ 10^* 
G. Indeed, it has been suggested that magnetic fields as strong as 10"*^^ G may be easily 
generated by dynamo processes in nascent neutron stars (Thompson & Duncan 1993), and 
may be required for soft gamma-ray repeaters (Paczynski 1992; Duncan &; Thompson 1993). 

Following the pioneering study of non-magnetic neutron star atmospheres by Romani 
(1987), recent atmosphere modeling has taken account of the transport of different photon 
modes through an ionized medium in strong magnetic field (Pavlov et al. 1994, 1995; Zavlin et 
al. 1995). Neutral atoms have been studied in detail for zero-field atmospheres (Rajagopal &; 
Romani 1996; Zavlin et al. 1996). Magnetic atoms have also been included in some models 
(Miller 1992), although many problems related to the treatment of the bound states in a 
strong magnetic field still remain. One outstanding issue concerns the non-trivial coupling 
between the center-of-mass motion and internal structure of the atom (e.g., Avron, Herbst 
&; Simon 1978; Herold, Ruder & Wunner 1981). This has been considered in detail in our 
recent paper (Lai & Salpeter 1995; hereafter Paper III) and used to derive the generalized 
Saha equation for ionization-recombination equilibrium. With this, we are now in a good 
position to calculate reliable composition and construct the complete phase diagram of the 
magnetic hydrogen surface of a neutron star. The published atmosphere models have shown 
that neutral hydrogen is important for zero-field, cool atmospheres (T ~ 10^ K), but we 
shall see that even at B12 ~ 0.01 its importance extends to higher temperatures. Moreover, 
we shall find that at high field strength polyatomic molecules (for B12 ^ 1) and even the 
condensed phase (for B12 ^ 100) are important. 

Our paper is organized as follows. In §2 we discuss the physics of various hydrogen bound 
states in a strong magnetic field. In §3 we study the physical properties of the nondegenerate 
atmosphere, including the abundance of various species. We consider the phase equilibrium 
of the metallic state in §4. Some astrophysical implications of our results are discussed in 
§5. Appendix B includes a study of the magnetic field induced pycnonuclear reactions that 
occur on neutron star surfaces. 

Throughout the paper we shall use real physical units and atomic units (a.u.) inter- 
changeably, whichever is more convenient. Recall that in atomic units, mass and length 
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are expressed in units of the electron mass nie and the Bohr radius = 0.529 x 10~^ cm, 
energy in units of 2 Rydberg = e^/oo = 2 x 13.6 eV; field strength in units of Bq (Eq. [1-1]), 
temperature in units of 3.15 x 10^ K, and pressure in units of e^/a^ = 2.94 x 10^^ dynes cm~^. 



2. HYDROGEN BOUND STATES IN STRONG MAGNETIC FIELD 

Here we briefly review the physics of various forms of hydrogen bound states in a strong 
magnetic field. The detailed quantum mechanical calculations are described in Paper I and 
Paper II, where extensive references can be found. Our main results are summarized in 
various fitting formulae for the binding energies and in Table 1. 

Our discussions are based on nonrelativistic quantum mechanics, even for extremely 
high magnetic field, B ^ B^ei = {hc/e^YBo = 4.414 x 10^'^ G, at which the transverse 
motion of the electron becomes relativistic. This nonrelativistic treatment of bound states 
is valid because (i) the electron remains nonrelativistic in the z-direction (along the field 
axis) as long as the binding energy is much less than mgC^, and (ii) the shape of the Landau 
wavefunction in the relativistic theory is the same as in the nonrelativistic theory, as the 
cyclotron radius 



1/2 



61/2 



(a.u.) = 2.57 x 10-^°bJ^^ (cm), 



is independent of the particle mass. Thus our results should be reliable even for the highest 
field strength considered in this paper. In §2.3, we shall include a discussion of the density- 
induced relativistic effect when we consider the metallic hydrogen state. 

2.1 Atoms 

In a superstrong magnetic field satisfying 6 >> 1, the spectra of the H atom can be 
specified by two quantum numbers (m, v), where m = 0, 1, 2, • • • measures the mean separa- 
tion pm = {2m + l)^^^p of the electron and proton in the transverse direction (perpendicular 
to the field), while is the number of nodes of the electron's z- wavefunction (along the field 
axis). The states with ly ^ resemble the zero-field hydrogen atom with small binding energy 
\E^\ ~ l/(2i/^) and we shall mostly focus on the tightly-bound states with u = 0. For the 
ground state (0, 0), the sizes (in atomic units) of the atomic wavefunction perpendicular and 
parallel to the field are of order L_i_ ~ /5 = b~^^^ and Lz ~ where I = \nb. The binding 
energy |i?(H)| (or the ionization energy Qi) of the atom is given by 



Qi = \E{a)\ ~ 0.16 P (a.u.) ~ 161 



ln(4265i2) 
In 426 



(eV). 



(2.1) 



The numerical factor 0.16 in equation (2.1) is an approximate value for B12 ^ 1 (More 
accurate fitting can be found in Paper II). Some numerical values are given in Table 1. The 
tightly-bound excited states (m, 0) have the transverse size L± ^ pm = [{2m + l)/b]^^^. For 
2m + 1 « b, we have ~ l/ln(l/p^), and 



Err, 0^ -0.16 It 



-0.16 In 



2m +1 



(a.u.) (for 2m + l «b). 



(2.2a) 
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For 2m + 1 ^ 6, we have 



pm , and the energy levels are approximated by 



0.6 




br 2m + 1^6). 



(2.26) 



Note that, unlike the field-free case, the excitation energy AE^ = \E(R)\ — \Ejn\ is small 
compared with |£'(H)|. 

The above results assume a fixed Coulomb potential produced by the proton (i.e., infinite 
proton mass). The use of a reduced electron mass memp/{me + nip) only introduces a very 
small correction to the energy (of order nig/nip). However, as mentioned in §1, the effect 
of the center-of-mass motion on the energy spectra is rather subtle in strong magnetic field. 
We will come back to this point in Sec. 3.2 (see Paper III). 



In a superstrong magnetic field, the mechanism of forming molecules is quite different 
from the zero-field case (Ruderman 1974; Papers I-II). The spins of the electrons of the atoms 
are all aligned anti-parallel to the magnetic field, and therefore two atoms in their ground 
states (m = 0) do not easily bind together according to the exclusion principle. Instead, 
one H atom has to be excited to the m = I state. The two H atoms, one in the ground 
state (to = 0) , another in the m = 1 state then form the ground state of the H2 molecule by 
covalent bonding. Since the "activation energy" for exciting an electron in the H atom from 
Landau orbital m to (to + 1) is small (see Eq. [2.2]), the resulting molecule is stable. In this 
way, more atoms can be added to form a larger molecule, in contrast to the field-free case. 

For a given magnetic field, as the number of H atoms N increases, the electrons occupy 
more and more Landau orbits (with to = 0, 1, 2,---,A'' — 1), but the length of the chain 
decreases so that the volume per electron of the electron distribution is of order (neglect- 
ing logarithmic factors) {bN)~^. Beyond some critical number Ng ~ [b/l'^y^^, it becomes 
energetically more favorable for the electrons to settle into the inner Landau orbitals (with 
smaller m) with nodes in their longitudinal wavefunctions (i.e., v 7^ 0). For N ^ Ns, the 
energy per atom in the Hjv molecule, \E(R]\f)\/N , asymptotes to a value 6^/^, independent 
of A'', and the volume per electron to ~ 6~^/^ (Paper I; see §2.3). For a typical magnetic 
field strength of interest here, the energy saturation point is A/'g ~ 3 — 5. 

The dimensions of the H2 molecule parallel and perpendicular to the magnetic field are 
comparable to those of the atom. The binding energies also approximately scale as (In 6)^, 
but they are numerically smaller than the ionization energy of H atom. More precisely, the 
dissociation energy of H2 can be fitted by 



= 2|£(H)| - |E(H2)| =0.106 [1 + r In (6/6crit)] (a.u.), r ~ O.U^'^, (2.3) 



with an accuracy of ^ 5% for 1 ^ B12 ^ 1000, where bcrit — 1.80 X 10* is to be defined in 
equation (3.8) (The superscript "(00)" implies that the zero-point energy of the molecule is 
not included in see §2.4). Thus Qi°°^ ~ 46 eV for 3^2 = 1 and Qi°°^ ~ 150 eV for 

B12 = 10 (cf. Table 1). By contrast, the zero-field dissociation energy of H2 is 4.75 eV. 

For the ground state of II2, the molecular axis and the magnetic field axis coincide, and 
the two electrons occupy the to = and to = 1 orbitals, i.e., (toi, TO2) = (0, 1). The molecule 
can have different types of excitation levels (Paper II): 

(i) Electronic excitations. The electrons occupy diflPerent orbitals other than (toi, TO2) = 
(0,1), giving rise to the electronic excitations. The energy difference between the excited 
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state (mi,m2) (with i/i = U2 = 0) and the ground state (0, 1) is approximately proportional 
to lnf>, as in the case for atoms. Typically, only the single-excitation levels (those with 
mi = and m2 > 1) are bound relative to two atoms in the ground states. Another type of 
electronic excitation is formed by two electrons in the (m, u) = (0, 0) and (0, 1) orbitals. The 
dissociation energy of this weakly-bound state is of order a Rydberg, and does not depend 
sensitively on the magnetic field strength. 

(ii) Aligned vibrational excitations. These result from the vibration of the protons along 
the magnetic field axis. For the electronic ground state, the energy quanta of small-amplitude 
oscillations is approximately given by 

nu\\ ~ 0.13 (Inbf/^iJ-^/^iai.u.) ~ 0.12 {\nbf/^ (eV), (2.4) 

where /i = mp/2me = 918 is the reduced mass of the two protons in units of the electron 
mass. Thus ^(jy ~ 10 eV at B12 = 1 and /iwy ~ 23 eV at B12 = 10, in contrast to the 
vibrational energy quanta ^Wvib — 0.52 eV for H2 molecule in zero magnetic field. 

(iii) Transverse vibrational excitations. The molecular axis can deviate from the mag- 
netic field direction, precesses and vibrates around the field line. Such an oscillation is the 
high-field analogy of the usual molecular rotation. The important difference is that here this 
"rotation" is constrained around the magnetic field line. For the ground electronic state, the 
excitation energy quanta hLO±o is given by 

nu±o ^ 0.125 6^/2 (In 6)//-^/2(a.u.) ~ 0.11 6^/^(ln6) (eV), (2.5) 

where the subscript "0" indicates that we are at the moment neglecting the magnetic forces 
on the protons which, in the absence of the Coulomb forces, lead to the cyclotron motion of 
the protons (see §2.4). Thus huj±o — 14 eV at B12 = 1 and huj±o = 65 eV at -B12 = 10. 

It is important to note that in strong magnetic field, the electronic and vibrational 
(aligned and transverse) excitations are all comparable, with huj±o ^ . This is in contrast 
to the zero-field case, where one has Aeeiec ^ ^^vib S> hoOj-ot- 

2.3 Infinite Chains and the Condensed Metallic State 

For N >> Ng. the binding energy per atom in a H^r molecule saturates, and the 
structure of the molecule is the same as that of an infinite chain Hoc. By placing a pile of 
parallel infinite chains close together (with spacing of order a.u.), a three-dimensional 

condensed metal can be formed (e.g., body-centered tetragonal lattice; Ruderman 1971). 
We offer our assessment on various calculations of the binding energy of magnetic Coulomb 
lattice. 

2.3.1 Uniform Electron Gas Model and its Extension 

The binding energy of the magnetic metal at zero pressure can be estimated using the 
uniform electron gas model (e.g., Kadomtscv 1970). Consider a Wigner-Seitz cell with radius 
rg, volume Vg = ATrr^/^ (the mean number density of electron is then = l/Vg). When 
the electron Fermi energy p|,/(2me) = {neh^c/2eB)^/2me is less than the cyclotron energy 
hUe, or when the density satisfies 

rie < riLandau = ^ J (2me?iu;e) = = 0-0716 63/2 (a.u.); 

he h V27r2p3 (2.6) 

'^s.Landau = 1-49 6"^/^ (^.U.), 
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the electrons only occupy the ground Landau level. In the simple uniform electron gas model, 
the energy per cell of the metallic hydrogen can be written as 

, , 37r2 0.9 , , , , 

^^(^^) = 85^6 (^-^O, (2.7) 

where the first term is the kinetic energy = (27r^/3me)(7i^c/eS)^ng, and the second 
term is the Coulomb energy. For the zero-pressure metal, we have dEg/drg = 0, and the 
equilibrium r<j and energy are then given by 

r,,o - 1.90 6-2/5 (a.u.), ^^,0 - -0.395 6^/5 (a.u.). (2.8) 

The corresponding mass density is p<j,o — 560 5^2^ g cm"'^, which is much smaller than the 
Landau density, defined in Eq. (2.6), by a factor 0.48 6"'^/^°. The pressure P in a metal 
compressed to, say, twice the zero-pressure density pgfi, is a few times r~Q\Es,o\, i.e., of order 

P ~ 0.1 6^/^ ~ 10^ 5^2^ (a.u.). Pressures in a neutron star surface layer of interest are much 
smaller than this and we need to consider only Ue < nLandau • 

We now discuss several corrections to the simple uniform electron gas model for metallic 
hydrogen. 

(i) Relativistic effect. As noted before, the use of non-relativistic equations for the 
transverse motion of the electrons is a good approximation even for B ^ i?rei — 137^.60. We 
can show that the density-induced relativistic effect is also small. The critical density nLandau 
for onset of Landau excitation is still given by equation (2.6). The relativistic parameter is 
Xe = pF/{mec) = (ne/nLandau)(2-B/Srei)^/2. At zcro-prcssurc density, using equation (2.8), 
we have Xe — 5.03 x 10~^b^^^ . Thus near zero-pressure density, the relativistic effect is always 
negligible for the range of field strength of interest in this paper. 

(ii) Coulomb exchange interaction. The exclusion principle for the electron results in an 
exchange correction to the Coulomb energy. The Dirac exchange energy (in atomic units) 
per electron is given by E^x = — 3F/(46rg), where F is a function of the ratio ne/riLandau 
(Fushiki, Gudmundsson & Pethick 1989). The effect of this (negative) exchange interaction 
is to increase rs,o and |£^s,o|. 

(iii) Non-uniformity of the electron gas. The Thomas-Fermi screening wavenumber Utf 
is given by (e.g., Ashcroft & Mermin 1976) kj,p = 4:iTe'^D{EF), where D{Ep) = due/dEp 
is the density of states per unit volume at the Fermi surface E = Ep = ]?|,/(2me). Since 
rig = {2eB/h'^c)pF, we have D{Ef) = ne/2EF = {mg/ne){2eB/h'^c)'^, and 

^TF=(^)'^%rf/2(a.u.). (2.9) 

(More details on the electron screening in strong magnetic field, including the anisotropic 
effect, can be found in Horing 1969). The gas is uniform when the screening length k^p is 
much longer than the particle spacing r^, i.e., kq^prg << 1. For the zero-pressure metal, 
using equation (2.8), we have kpprg ~ 1.83, independent of B. Thus even as J5 ^ 00, 
the electron non-uniformity must be considered for the zero-pressure metal. This effect can 
be studied using the Thomas-Fermi type statistical model, including the exchange and the 
Weizsacker gradient corrections (see Fushiki et al. 1992 and references therein). 

2.3.2 Cohesive Energy of the Condensed State 

Although the simple uniform electron gas model and its Thomas-Fermi type extensions 
yield reasonable binding energy for the metallic state, their validity and accuracy in the strong 
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field regime cannot be easily justified, and it is the difference between the energy of the metal 
and that of the atom that determines whether the metal is bound. One uncertainty concerns 
the lattice structure of the metal, since the Madelung energy can be very different from the 
Wigner-Seitz value (the second term in Eq. [2.7]) for a non-cubic lattice. In principle, a three- 
dimensional electronic band structure calculation is needed to resolve this problem, as Jones 
(1985, 1986) has attempted for carbon and iron using density functional theory. However, 
the electron correlation has not been taken into account, and, like the Thomas-Fermi type 
statistical models, the accuracy of density functional approximation in a strong magnetic 
field is not yet known (for a review of density functional theory as applied to non-magnetic 
terrestrial solids, see, e.g., Callaway & March 1984). 

We consider the self-consistent Hartree-Fock method (Neuhauser, Koonin & Langanke 
1987; Paper I) to be the most reliable approach to the problem. This has only been done for 
one-dimensional chains. Our numerical results for the energy (per atom) of the Hoo chain 
can be fitted (to within 2% accuracy for B12 up to ~ 10"^) to a form similar to equation (2.8): 

= -0.76 6°-^^ (a.u.) = -194 B%^^ (eV). (2.10) 

Note that lE'ool in equation (2.10) is larger than equation (2.8) by a factor of 1.8. The 
cohesive energy of the Id chain (energy release in H-I-Hoo =Hoo_|_i) is given (to within ^ 10% 
accuracy) by 

Q(^) = |£oo|-|^(H)|~0.76 6°-^^-0.16(ln6)2 (a.u.) (2.11) 

where the superscript "(00)" indicates that the zero-point energy correction to the cohesive 
energy has yet to be included (§2.4). 

The energy difference AEg = |-Es,o| — l-E^ool between the 3d metal and the Id chain 
must be positive and can be estimated by considering the interaction (mainly quadrupole - 
quadrupole) between the chains. We have found (Appendix A) that this difference is probably 
between 0.4% and 1% of l-Eool- Therefore, for hydrogen, the binding of the three dimensional 
metal results mainly from the covalent bond along the magnetic field axis, not from the chain- 
chain interaction. The cohesive energy of the zero-pressure H metal is Qg = Qoo + . 

A note about the difference between iron and hydrogen is appropriate at this point. For 
Fe, it is has been found that at B12 ~ 1 — 10, the infinite chain is not bound relative to the 
atom (Jones 1985; Neuhauser et al. 1987), contrary to what the original calculations (e.g.. 
Flowers et al. 1977) indicate. Therefore, the chain-chain interaction must play a crucial role 
in determining whether the three dimensional zero-pressure Fe metal is bound or not. The 
main difference between Fe and H is that for the Fe atom at B12 1, many electrons are 
populated in the 7^ 1 states, whereas for the H atom, as long as 6 >> 1, the electron always 
settles down in the = tightly-bound state. Therefore, the covalent bonding mechanism 
for forming molecules (cf. §2.2) is not effective for Fe at B12 ~ 1- However, for sufficiently 
high B field, when Uo/Z » v^2ZTTp, or B12 » 100(.^/26)^, we expect the Fe chain to 
be bound in a similar fashion as the H chain discussed here. 

2.4 Zero- Point Energies and Relative Binding Energies 
of Different Forms of Hydrogen 

For the H2 molecule, the zero-point energy of the aligned vibration is huj\\/2 (Eq. [2.4]). 
Equation (2.5) for the zero-point energy of the transverse oscillation includes only the con- 
tribution of the electronic restoring potential ij,u;']_qR']_/2. Since the magnetic forces on the 
protons also induce a "magnetic restoring potential" ^u)pB?^/2, the zero-point energy of the 
transverse oscillation is (Paper 11) 

nLO± = hiuio + uD^/^ - hcop, (2.12) 
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where hiOp = TieB/{mpc) ~ 6.3 eV is the cyclotron energy of proton. Thus dissociation 
of H2 taking into account the zero-point energies is 

Q2 = Q^r^ - {]^nuj\\ + nuj^ . (2.13) 

Now consider the zero-point energy in the metallic hydrogen. Neglecting the magnetic 
force, the zero-point energy Ezp of a proton in the lattice is of order HQp, where Qp = 
(47re^ne/mp)^/^ is the ion plasma frequency. Thus for the zero-pressure metal, we have 

E^p ~ hQp = 0.040 r7^/^ (a.u.) ~ 0.015 6^/^ (a.u.). (2.14) 

This is much smaller than the total binding energy (Eq. [2.10]) unless B12 ^ 10^. This means 
that for the range of field strength considered in this paper, the zero-point amplitude is small 
compared to the lattice spacing. Thus quantum melting is not effective (Ceperley & Alder 
1980) and the metal is a solid at zero temperature. Accurate determination of E^p requires a 
detailed understanding of the lattice phonon spectra. At zero-field, Monte-Carlo simulations 
give Ezp — 3h^}prj/2, with rj ~ 0.5 (Hansen & Pollock 1973). For definiteness, we will adopt 
the same value for Ezp in a strong magnetic field. Taking into account the magnetic effect 
on the proton, the corrected cohesive energy of Id chain is expected to be 



(00) 
00 



1 ( uj'^ 

-hQpTj + ni^ + rjHl 




(2.15) 



with T] ~ 0.5. As mentioned in §2.3, the relative binding energy between 3d condensate 
and Id chain is close to 0.4% — 1% of \Eoo\- We shall consider the condensed phase only 
for B12 > 10, where the ratio Qoo/|-E^oo| ranges from 0.3 to 0.7 for B12 = 10 — 500. For 
definiteness we shall express the cohesive energy of the 3d condensate as 

Q, = \E,,o\ - \E{R)\ = Qoo + AE, = (1 + a)Qoo, (2.16) 

with a ~ 0.01 - 0.02. 

To summarize, we plot in Figure 1 the energy releases Qi, Q2 and Qoo for e-|-p=H, 
H-|-H=H2, and H+Hqo =Hoo+i respectively. Some numerical values are given in Table 1. 
The zero-point energy corrections for Q2 and Qoo have been included in the figure (if they 
are neglected, the curves are qualitatively similar, although the exact values of the energies 
are somewhat changed.) Although b » 1 satisfies the nominal requirement for the "strong 
field" regime, a more realistic expansion parameter for the stability of the condensed state 
over atoms and molecules (cf. Eqs. [2.1], [2.3], [2.8] and [2.11]) is the ratio fe°-'^/(ln6)2. This 
ratio exceeds 0.3, and increases rapidly with increasing field strength only for b ^ 10^. For 
B12 ^ 10, we see that Qi > Q2 > Qoo (cf. Fig. 1), the condensed metallic phase may (or 
may not) exist in the atmosphere. Fortunately, at the temperatures of interest in this paper 
we shall need to consider the condensed state only for B12 ^ 10, where Qi > Qoo > Q2 (for 
10 <, B12 ^ 100) or Qoo > Qi > Q2 (for B12 ^ 100). This will have important consequences 
on the composition of the saturated vapor above the condensed phase (§4). 

2.5 Surface Energy of Droplets and Large Molecules 

For a single linear molecule Hjv we write Qn as the dissociation energy into Hjv-i+H 
and denote the limit for an infinite chain as Qoo- Numerical values and fitting formulae for 
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TABLE 1 



Relative Binding Energies (in eV) of H, H2 and H, 







q1°°^ 


^ z 


^ 00 


0^ 


0.1 


76.4 


14 


13 






0.5 


130 


31 


21 






1 


161 


46 


32 


29 


20 


5 


257 


109 


80 


91 


71 


10 


310 


150 


110 


141 


113 


50 


460 


294 


236 


366 


306 


100 


541 


378 


311 


520 


435 


500 


763 


615 


523 


1157 


964 


1000 


871 


740 


634 


1630 


1350 



NOTE: B12 = -6/(10^^ G). Qi is the ionization energy of H atom, Q2 is the disso- 
ciation energy of H2, and Qoo is the cohesive energy of Id H chain. The superscript 
"(00)" implies infinite proton mass. The given numerical results are generally accu- 
rate to within 10%. Calculations of Qoo are not reliable for B12 ^ 1. 



Q2 and Qoo are presented in §§2.2-2.4. Results (without zero-point energy correction) for 
A'^ = 3, 4 have also been given in Paper I. While accurate energies for larger molecules are 
needed in order to establish the proper scaling of Qn as a function of and B, we note the 
following trend: The ratio Q2/Q00 (and to some extent Q3/Q00) decreases appreciably with 
increasing B, but Qn/Qoo varies little with B for large A^. For a given B, the ratio Qn /Qoo 
reaches a maximum near N ^ Ng (the saturation point) and then approaches unity from 
above for >> A^^g. 

For the phase equilibrium between the condensed metal and Hjv molecules in the vapor 

(§4) we will need the "surface energy" Sn, defined as the energy release in converting the 
3d condensate H^^oo and a H^v molecule into Hg^oo+Ar- Clearly Si = Qg = (a + l)Qoo is 
the cohesive energy defined in §§2.3-2.4. For linear Hjv molecule with energy (per atom) 
En = -E(Hjv), we have 

Sn = NiEN - Eg) = NAEg + N{En - ^co) = NaQ^o + fQoo, (2.17) 

where the first term on the right-hand side comes from cohesive binding between chains, and 
the second term is the "end energy" of a Id chain. For N > 2, we can express Sn in terms 
of Qjv's via 

N r N ' 

Sn = NQg-Y,Qi = NaQoo+ Qoo + Y.^Qoo - Qi) 

The dimensionless factor / in Eq. (2.17) is of order unity: For AT = 4, it equals 0.23, 0.6, 1.1 
and 1.6 for B12 = 1, 10, 100 and 500, respectively. The asymptotic values of / for A'' >> A^^ 
should be close to these values (since Qn ^ Qoo ioi N >> Ng). 

Since a ~ 0.01, representing the cohesion between infinite chains, is so small, the first 
term in equation (2.17) is appreciable only when N ^ ~ 100. In that case, the molecular 



(2.18) 
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configuration which minimizes the surface energy Sn of the condensed metal is not a linear 
chain, but some highly elongated "cylindrical droplet" with A^_l parallel chains each contain- 
ing A''|| atoms with iVj_A''|| = A''. For such a droplet, the "end energy" is of order A^_l/Qoo- On 
the other hand there are ~ VN± "unpaired" chains in such a droplet, each gives an energy 
N\\aQoo- Thus the total surface energy is of order [A^_l/ + a{N/N±)y/N±] Q^o- The mini- 
mum surface energy Sn, for a fixed N ^ 2fa~^ = Nc, is then obtained for A^_l ~ (N/Nc)"^^^, 

/ ^ \ 2/3 2 f 

SNC^3fl—j Qoo, for AT ^ iV, = ~ 200. (2.19) 

Thus, although the optimal droplets are highly elongated, the surface energy of the condensed 
metal still grows as {N/200f/^ for N ^ 200. 

2.6 Other Species: H" and 

For completeness, here we also give fitting formulae for the binding energies of the 
negative ion H~ and the molecular ion Hg". The energy release in H-|-e=H~ is 

g(H-) ~ 0.014 (In 6)2 (a.u.) (2.20) 

Thus the ionization energy of H~ is 13 eV for B12 = 1 and 24 cV for B12 = 10, as compared 
to its value of 0.75 eV for B = 0. The energy release in H-|-p=H^ is 

Q(H+) ~ (0.008 + 0.11 \nb){\nbf (a.u.) (2.21) 

Equations (2.20)-(2.21) are both accurate to within 20%. 

Because of the small binding energies, both H~ and H2 are likely to have negligible 
abundance in a typical neutron star atmosphere. However, the H~ ion might, in principle, 
contribute appreciably to the atmospheric opacity. 



3. THE SAHA EQUILIBRIUM AND WARM ATMOSPHERES 

We now consider the physical conditions and chemical equilibrium in neutron star atmo- 
spheres with photospheric temperature in the range Tph ~ 10^ — 10^ K and magnetic field 
strength in the range B12 ~' 0.1 — 20. These conditions are likely to be satisfied by most ob- 
servable neutron stars. For such relatively low field strengths, the atmosphere consists mainly 
of ionized hydrogen, H atoms and small H^v molecules, and we can neglect the condensed 
metallic phase in the photosphere. Although the density scale height of the atmosphere is 
only h ~ hTpii/iripg ~ 0.08Tph,5 (7f/ (cm), where Tph, 5 = Tph/(10^ K), g is surface gravita- 
tional acceleration, gi4 = (7/(10^'* cm s"^), the atmosphere has significant optical depth. In 
§4 we shall consider more extreme situation (S12 >> 10) when the non-degenerate atmo- 
sphere has negligible optical depth and the condensed metallic phase becomes important. 

Although the surface layers we consider can lie somewhat below the photosphere, we 
shall restrict to densities smaller than the zero-pressure condensate density (cf. §2.3) 

p«,o^ 560Sf/' (gcm-3). (3.1) 
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Thus the density is always below the critical "Landau density" (cf. Eq. [2.6]) 

PLandau = 7.08 X 10^ /^.B^i^ (g cm-^), (3.2) 

where He is the mean molecular weight per electron (/ig = 1 for ionized hydrogen). For our 
ranges of B and Tph we also have the inequality T « TiLOe/k = 1.34 x 10^i?i2 K, so that all 
electrons settle in the ground Landau state. The electron Fermi temperature Tp = Ep/k is 
then given by 

TF = 2mB^i{j-^ (K), (3.3) 
where p is the mass density in g cm~^. At the density of the condensate (Eq. [3.1]), we would 

^ 2/5 

have Tp ~ 8.4 x 1Q^B^2 K; but for p much smaller than ps^Oj we have Tph >> Tp, and the 
electrons are nondegenerate.^ 

3.1 The Photospheric Density and Pressure 

In a star with an interior energy source (such as an isolated cooling neutron star), the 
photon flux F = (JsbT'eff constant to very far below the photosphere (o"sb is the Stefan- 
Boltzmann constant). With an accretion energy source, this constancy holds only above 
the stopping layer for the infalling matter. However, the optical depth Tstop in this layer is 
much larger than unity for a warm atmosphere^, and we shall consider only optical depth 
r less than Tstop) where r = J ndy with k the Rosscland mean opacity and y = j pdz the 
column density. With the flux F = ashT^s assumed constant, a simple approximation for 
the temperature profile T{t) is 

r^(r)c.^r,1,0+r), (3.4) 

and the photospheric temperature Tph = Teff- (with Tph chosen as 2/3). With the gravitational 
acceleration g in the atmosphere constant, hydrostatic equilibrium gives for the pressure 
P = gy and, rewriting t = Ry, we have 

P(t) = |t; Pph:.|. (3.5) 

Note that Eq. (3.5) holds irrespective of the equation of state, even when there is a phase 
transition at some pressure (see §4). For the warm atmospheres considered in this section, 
the ideal gas law applies and the photospheric density is given by Pph — Pph^^p/ (kTph). 

For the temperature and density relevant to the atmosphere, free-free absorption dom- 
inates over electron scattering. In the zero-field case, assuming Kramer's opacity k(0) = 
KopT~^'^, with Ko c± 7.4 x 10^^ (cgs), we would obtain the photosphere density-temperature 

^ Equation (3.3) should be compared with the Fermi temperature in field-free case, Tp{0) = 
3.0 X 10^{p/pe)'^/^ (K). Clearly, high magnetic field lifts the degeneracy of electrons even at 
relatively high density. 

^ The column density i/stop of the stopping layer typically corresponds to a few to tens 
of Thompson depths (depending on the field strength), i.e., ystop ~ lO/Kesi where k^s is 
the electron scattering opacity (cf. Nelson et. al. 1993). Thus for an absorption-dominated 
atmosphere, Tgtop ~ 10Kfr//«es >> 1- 
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relation Pph(O) ~ 7 x 10~^ gl{^ '^ph\ gcm"'^. In a strong magnetic field, the radiative opacity 
becomes anisotropic and depends on polarization (Camito ct al. 1971; Lodenquai et al. 1974; 
Pavlov et al. 1994). For photons with polarization vector perpendicular to the magnetic 
field (the "extraordinary mode") the free-free absorption and electron scattering opacities 
are reduced below their zero-field values by a factor of (w/wg)^, while for photons polarized 
along the magnetic field (the "ordinary mode"), the opacities are not affected. Silant'cv and 
Yakovlev (1980) have calculated the appropriate average Rosseland mean free-free opacity. 
In the field and temperature regime of interest, an approximate fitting formula is 



k(S)~ 400/3 — k(0), (3.6) 

\hljJe ) 

with /3 c± 1. The resulting photosphere density and pressure are given by 

Pph(5)^0.5/rV2<^if TpYji^is (gcm-^), 
Pph(i?) ^ 4 X 10^2 /J-V^^jfTp^/^^ (dynecm-2). 

In Figure 2, this photosphere (r ~ 2/3) condition is shown in a temperature-density diagram. 
For /3 ~ 1/400, equation (3.7) also approximately characterizes the physical conditions of 
the deeper layer where the extraordinary photons are emitted. This "r^ ~ 2/3" line is also 
shown in Figure 2. Solving Eqs. (3.4)-(3.5) (more precisely, dP/dr = g/n), we can obtain the 
temperature profile of the atmosphere as a function of density. Some of such T — p profiles 
for different values of Teff = Tph are depicted in Figure 2. Clearly, T rph/2^/^ = 0.84rph 
as r — > 0, while T oc t^^^ oc p'*/^ as r ^ oo. 

Other sources of opacity such as bound-free and bound-bound absorptions will increase 
the opacity and reduce the photosphere density, but the above estimates define the general 
range of the physical parameters in the atmosphere if Teg is large enough for the neutral 
H abundance to be small. This rough estimates agree reasonably well with more detailed 
calculations of Pavlov et al. (1995). 

3.2. Ionization Saha Equilibrium 

We now consider the ionization-recombination equilibrium of the H atom. Previous 
treatments of this problem (e.g., Khersonskii 1987; Miller 1992) have assumed that the H 
atom can move across the magnetic field freely. This is generally not valid for the strong field 

regime of interest here. A free electron confined to the ground Landau level, the usual case 
for 6 >> 1, does not move perpendicular to the magnetic field. Such motion is necessarily 
accompanied by Landau excitations. When the electron combines with a proton, the mobility 
of the neutral atom across the field depends on the ratio of the atomic excitation energy 
(~ In 6) to the Landau excitation energies hcOp = heB / {ijipc) for the proton. It is convenient 
to define a critical field strength .Bcrit via 

6crit = ^ In 6cHt = 1-80 X 10^ Sent = ^crit^o = 4.23 X 10^3 G. (3.8) 

Thus for B ^ i?crit ; the deviation from the free center-of-mass motion of the atom is signifi- 
cant (Paper III). 

In a strong magnetic field, the center-of-mass motion of the atom is specified by a pseu- 
domomentum K, whose component along the field axis is simply the usual liner momentum. 
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while the perpendicular component K_l measures the mean transverse separation of the elec- 
tron and proton (Avron, Herbst & Simon 1978; Herold, Ruder & Wunner 1981; Paper III 
and references therein). For a H atom in the ground state, the total energy (in atomic units) 
is approximately given by (Paper III) 



(3.9) 



where we have defined 



M,.M(l + ^).M(l + il), 1^,^.0.64^6(1 + ^ 



(3.10) 



and M = mp + me ~ mp/m^ (a.u.), ^ ~ 2.8, |-E(H)| ~ 0.16P (a.u.) (cf. Eqs. [2.1]-[2.2]). 
Note that for « Kc, the dependence of So on K_i becomes Kj_/{2M±). Thus M± 
represents the effective mass for the transverse motion across the magnetic field. Potekhin 
(1994) has given more accurate numerical results for a few selected field strength, but the 
approximation in equation (3.9) is adequate for our purpose. 

Let the proton (free or bound) number density in the gas be Ug. The partition function 
(in atomic units) for a H atom in a volume Vg = AirVg/S = l/n^ is then 



Z(H) :^ Vg 



fMT 



1/2 



(3.11) 



where 



Mi 



_ 2M_LrY 



(3.12) 



for T << K'^/{2M±), and M'j_/M± approaches a constant less than unity when T becomes 
comparable to K^/{2M±). The partition function .Z'(H) has the same form as the zero-field 
expression except for the factor M'j_/M. In equation (3.11), z(R) = Zm(H)z,^(H), where 
Zm(ii) is the internal partition function associated with the m > excited states: 



.(H) ^ ( 



1 + e- 



-b/MT'^ 



"'-max T,r/ 

E^exp 



m=0 



M' 



(o.l6l^ + Ern + m^ 



(3.13) 



where M'j_^ is of the same order of magnitude as Af^, the factor (l + g"''/^^^) comes from 
the proton spin effects, and mmax is set by the condition L± ^ Vg or ^ Vg (cf. §2.1). The 
internal partition function Zi, associated with the u > states is close to unity (see Paper 
III for detail). 

The partition functions for the ionized states (free electron and proton) can be eas- 
ily obtained. With the atomic bound-state partition function given above, we can obtain 
the generalized Saha equation for the ionization-recombination equilibrium in the pres- 
ence of strong magnetic field. In the density and temperature regimes of interest, with 
T ^ K^/{2M±) ~ 0.32 / (1 + Ml/^b), we have 



n(H) 



b_ 
2^ 



M': 



J. N 1/2 

— I tanh 
27r J 



2MT 



tanh 



2T 



exp 



Qi 

T 



z{R), 



(3.14) 
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where n(H), Up and ng are the number densities of the different species. 



3.3 Dissociation Equilibrium of Molecules 

Accurate treatment of the dissociation equilibrium for Hjv molecules is complicated. 
However, since the molecular excitation energies are comparable to the excitations in H 
atom (§2.2), we expect that when T ^ K^/{2Mj_) ~ / (1 + Ml/^b), we can similarly use an 
effective mass description for the motion of the molecule across the magnetic field. As an 
estimate, we assume that the effective mass of Hjv molecule is M_L(Hiv) ~ A^M(1 + ^h/Ml). 
We similarly introduce the correction factor M'j_/M± given in equation (3.12). The Saha 
equation for the equilibrium process Hjv+H=Hjv+i + Qn+i is then given by 

n(HAf+i) f N + iV^'^ f MT\~^^^ / M \ fQN+i\ z(R 



n(Hjv)n(H) 



fN + iy^yMTy^'f M\ [Qn+a ^(h^+i) 



where z(H]\f-^-i), z(Hjv), ^(H) are the internal partition functions of Htv+i, H^t, H respec- 
tively. To estimate the internal partition functions of the molecules, one need to include 
various molecular excitation levels as discussed in §2.2. But presumably the ratio of these 
internal partition functions is of order unity (this approximation is adopted in the calcula- 
tions presented in §3.5), since the molecular excitation level spacing is not necessarily smaller 
than the electronic excitation (§2.2). 

3.4 Non-ideal Gas Effect 

So far we have assumed the gas to be sufficiently dilute to be treated as an ideal gas. 
When the gas density increases, the interaction between particles becomes important. The 
effect of the atom-atom interaction potential Ui2{r) on the ionization equilibrium is to modiiy 
the chemical potential of the atomic gas by the amount 

A/i(H) = n(H)A;r J d^r{l - e^'^/''^), (3.16) 

(e.g.. Landau & Lifshitz 1980). For r much larger than the atomic length ~ 1/Z of the 
elongated atom, the potential U12 is due to the quadrupole-quadrupole interaction U12 ~ 
Q^(3 — 30 cos^ 0-1-35 cos^ d)/r^, where 9 is the angle between the vector r and the z-axis, and 
Q ~ eL^ is the quadrupole moment of the atom. Since the integration over the solid angle 
J dilUi2 = at large r, the contribution to A//(H) from large r is negligible. The main effect 
is then the "excluded volume effect": Let Va be the volume of the (highly non-spherical) 
electron distribution for one atom, so that we can set U12 00 due to quantum mechanical 
repulsion when two atoms overlap. We then have A^(H) ~ kTn{}l)va- The atom-atom 
interaction introduces a factor exp(— A///fcT) ~ exp(—nva) to the right-hand-side of the 
Saha equation (3.14). Therefore when the mean density in the gas becomes comparable to 
the internal density of the atom, the atom is mostly ionized — so called pressure ionization. 
One can similarly consider the interaction between a charged particle (electron or proton) 
and the atom, given by the quadrupole-monopole potential U12 ~ eQ{3 cos^ 9 — 1) /r^ . Again 
since / dQUi2 = 0, the long-range term is negligible. We thus obtain a similar factor of order 
^ exp{—nVa) correction to the Saha equation. 

The "atomic volume" Va is of order ttL^L^, where L± ~ is the radius perpen- 

dicular to the z-axis and ~ l/l the length in the z-direction. The ratio of this volume 
to the volume per electron in the condensed state, 47rrg o/3 with Vg^o given in equation 
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(2.8), is then ~ 0.1 b^^^ /Inb, which increases slowly with increasing b and is of order 0.1 for 
Bi2 ~ 1 — 10'^. The "excluded volume" of the atom may be larger due to its elongated shape, 
but Va ^ Ll ^ 47rr^ Q/3. Thus we have A//(H) ~ kT (p/ps^), where ps,o given by Eq. (3.1). 
The corresponding pressure-ionization factor in the Saha equations is exp(— A/^/fcT). 

3.5 Results 

To get a qualitative overview of the relative abundances of different forms of hydrogen, 
we have determined the H atom half-ionization and H2 half-dissociation curves in the T — p 
diagram. In Figure 2, we show the results for B12 = 1 and B12 = 10; we also plot the electron 

degeneracy line (T = Tp) and the typical photosphere conditions (§3.1). The half-ionization 
line is determined by setting rie/rig = Up/ug = n(H)/ng = 1/2, while the half-dissociation 
line is obtained with n(H)/ng = 1/2 and n(H2)/ng = 1/4 (so that H and H2 have equal 
mass fractions). Note that the gas density is simply p = mpUg, and the pressure ionization 
factor discussed in §3.4 is not included in the calculations. To the left of the half-ionization 
line (labeled "e-|-p=H") the atmosphere is mostly ionized, while to the right of the half- 
dissociation line (labeled "H-|-H=H2"), most of the material is in the molecular states, H2, 
H3, etc. Between the ionization line and the H2 dissociation line, the dominant species is H 
atom. From Eqs. (3.7) and (3.14), we can obtain the photospheric "ionization temperature", 
^phjirj where the ionization-recombination equilibrium gives 50% H and 50% e+p in the 
photosphere (r = 2/3). For B12 ~ 1 — 50, an approximate fitting formula is 

logio Tph.ir ^ 5.4 + 0.3 logio B12. (3.17) 

If the actual photospheric temperature Tph is greater than Tph,ir, then there is little neutral 
H near and above the photosphere (r ^ 1), and the neutral abundance increases slowly 
with depth. If Tph is only slightly smaller than rph,ir, then neutral H dominates in layers 
near the photosphere (r ~ 1), but (e+p) dominates at higher levels (r << 1), because T is 
almost constant for r ^ 1 while p decreases to zero as r ^ (cf. Fig. 2). For B12 = 1, the 
abundances of molecular species are negligible unless the photospheric temperature drops 
well below 10^ K. For B12 = 10, there exists a large amount of H2 in the photosphere when 

The "neutrality fraction" n(H)/ng as a function of temperature at a fixed density p = 
ngTUp = 0.01 g/cm^ and 0.1 g/cm^ is shown in Fig. 3. Results for a wide range of field 
strengths {B12 = 0.01, 0.1, 1, 10) are given and compared with the zero-field limit. In 
general, the neutral fraction is not a monotonic function of B, because it is determined by 
two opposite effects [cf. Eq. (3.14)]: the ionization energy Qi increases with increasing B, 
this tends to increase n(H)/ng; on the other hand, the phase space of free electron and 
proton, proportional to 6^/ tanh(6/2MT), also increases with increasing B, and this tends 
to decrease n(H)/ng. It is also easy to understand from Eq. (3.14) that the neutral fraction 
is in general a non-monotonic function of T at a given field strength. As seen from Fig. 3, 
even a relatively small field increases the neutrality appreciably above the field-free value, 
partly because Qi is larger, but also because the internal partition function Zm{H) is larger. 
Thus it may well be that & B = 10^ — 10^° G atmosphere (which characterizes millisecond 
pulsars) is very different from a field-free atmosphere (Rajagopal Sz Romani 1996; Zavlin et 
al. 1996). 

Figure 4 shows the fraction (by number) of various hydrogen gas species as a function 
of Tph at a point with P = O.SPph (cf. Eq. [3.7] with /3 = 1 and 514 = 1), r = 1/3 
and T = (3/4)i/'*Tph = 0.93Tph. Three magnetic field strengths {B12 = 1, 10, 100) are 
considered. The total number density ntot of free particles in the gas, i.e., ntot = Ue + Up + 



16 



n(H) + n(H2) + • • •, is given by ntot = P/kT = 1.6 x lO^^BisTpY 5 (cm"^). Neglecting the 

trace ion species such as H~ and Hj" (and the condensed phase), the abundances satisfy the 

condition 2Xp + X(H) + ^(Ha) + ^(Hg) H = 1, where Xp = Up/ntot and X(Hjv) = 

n(HAr)/ntot- For determining Xp, X(R), X{R2), the values of Qi, Q2 given in Table 1 
are used, while for determining the fractions of larger molecules, the following approximate 
ansatz is adopted in our calculation (cf. §2.5): For H3 and H4, the ratio Q^^^ /Q^^ given in 
Paper 1 is assumed, and the zero-point energy correction to Qn is introduced by setting = 
(Qiv°VQoc°'')Qoo with given in Table 1; For A'' > 4, where no numerical result is currently 
available, we simply set Qjv = Qoo- Clearly, the uncertainty is largest for intermediate-sized 
molecules (N Ns), but we have checked that moderate variations (less than 30%) of 
these intermediate Qjv's do not change our results significantly. The approximation should 
become better for N » Ng as Qn asymptotes to Qoo- For B12 = 1, 10, we see that the 
atmosphere is dominated by atoms and relatively small molecules, and the total gas density 
p = mpUg = ?Ti,p[np + n(H) +2n(H2) +3n(H3) + • • •] remains much less than the condensation 
density Ps,o- For the extreme field strength B12 = 100, which we include here for illustrative 
purpose, increasingly large molecules dominate the atmosphere as Tph decreases, and the gas 
density approaches p^^o- This is indicative of a phase transition which will be discussed in 
§4. 

For layers deeper than the photosphere, i.e., for optical depth r increasing to large values, 
the ionization fraction changes slowly since the temperature-density track in the atmosphere, 
p oc T^/^ (for r >> 1), is close to the track of constant (e+p)/H ratio. Indeed, from Fig. 2 
we see that the ionization fraction varies by only a factor of a few below the photosphere. 
The zero-pressure density psfi of the condensed metal is given by Eq. (3.1) and shown by 
the heavy vertical lines in Figure 2. Note that the half-ionization lines and half-dissociation 
lines arc meaningful only in the density regime to the left of the condensation density lines 
and the degeneracy lines. When r is sufficiently large that p approaches Ps^, non-ideal 
gas efi!'ects and/or condensation set in. With T oc p^^^, the temperature at p ~ ps,o is of 
order T ~ 20 Tph. For a warm atmosphere this exceeds the critical temperature for phase 
transition, Tcrit ^ 10^ K (see §4). Thus there is no distinct condensed phase, although 
the atoms would overlap at p ^ ps^ and be pressure-ionized. The electrons form a uniform 
fluid, but stay in the ground Landau level until p exceeds pLandau (the dot-dashed vertical 
lines in Fig. 2). 



4. ULTRAHIGH FIELDS AND COOL ATMOSPHERES: 
THE CONDENSED PHASE 

As discussed in §2, the cohesive energy Qg and surface energy Sn of the condensed 
metallic hydrogen increase more rapidly with increasing fleld strength than the ionization 
energy Qi of H atom and the dissociation energy Qn of small Hjv molecules (cf. Fig. 1). We 
shall see below (cf. Eqs. [4.5]-[4.7]) that Qg = Si, Sn and {Qg + Qi)/2 determine the partial 
saturation vapor densities of H, and (p+e) respectively. For B12 » 10, the cohesive 
energy is greater than Qi and Q2, and is close to \Eg^Q\ ~ 8B^^ a.u. In this case there may 
exist a critical temperature Tcrit, below which there is a first order phase transition between 
the condensed metallic hydrogen^ and the gaseous vapor. We explore the possibility of this 

^ The condensed metal may be solid or liquid. Using the zero-field criterion as an estimate, 
the melting temperature is determined by F = e^/^rgkT) = IS.TB^^^ /T^ c± 180, where we 
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phase separation and the partial saturation vapor pressure of various constituents in this 
section. 



4.1 Phase Equilibrium Conditions 

Consider the equilibrium between the condensed metallic phase (labeled by the subscript 

"s") and the non-degenerate gaseous phase (labeled by "g"). The electron number density 
in the metallic phase is Ug = l/Vs = 3/(47rrg), where is the mean proton spacing (the 
radius of Wigner-Seitz cell). The gaseous phase consists of a mixture of free electrons, 
protons, bound atoms and molecules, with the total baryon number density Ug = 1/Vg = 
np + n{}i) +2n(K2) + ■ ■ ■■ Wc neglect the small concentration of H~ , and other molecular 
ions in the gas, thus we have rie = Up. Phase equilibrium requires the temperature, pressure 
and the chemical potentials of different species to satisfy the conditions: 

T,=Tg= T, 

Ps = Pg = [2np + n(H) + ^(Ha) + n(H3) + • • ■]kT = P, ^^-^^ 
Us = yu(H) = n^ + i_ip = -//(Ha) = -//(H3) = • • • 

First consider the equilibrium between the condensed metal and H atoms in the gaseous 
phase. Near the zero-pressure metal density, the electron Fermi temperature is Tp.o — 
8.4 X 10^3^2^ K (cf. Eqs. [3.1], [3.3]). Thus at a given temperature, the metal becomes more 
degenerate as B increases. Let the energy per Wigner-Seitz cell be Es{rs)- The pressure and 
chemical potential of the condensed phase are given by 

1 dE 

Ps = --r^-r^, lis = E,{r,)+P,V,c^E,^o + PsVs,o, (4.2) 
47rrj avs 

where the subscript "0" indicates the zero-pressure values, and we have assumed that the 
vapor pressure is sufficiently small so that the deviation from the zero-pressure state of 
the metal is small, i.e., 5 = [(r^ — r<j^o)/''s,o| << 1- This is justified when the saturation 
vapor pressure Pgat is much less than the critical pressure Pcrit for phase separation, or 
when T « Tent- The finite temperature correction A/Xg to the chemical potential of the 
condensed phase is given by ^ 

- ^ 0.10Br2'/^r5, A;..(r) 0.85Sra^/^r| (eV), (4.3) 

where the Fermi energy Ep = kTp = 97r^/(86^rf) ~ 0.236 6^/^ (a.u.), and for the density we 
have used the zero-pressure value. The partition function of H atoms is given by equation 
(3.11). Using n(H) = exp[//(H)/T]Z(H)/T^ and the equilibrium condition /x(H) = /Xg, we 
obtain the number density of H atoms in the saturated vapor: 



.MTV'^(M'^\ [~Q,, + PV:,.o + A/y„(r)-A//(H) 



(4.4) 



5 /2 

have used the zero-pressure value for given by Eq. (2.8). Thus for B12 ^ 287 the 
zero-pressure metal is likely to be a liquid. 

For kT « huJe and Ep << huje so that electrons arc all in the ground Landau level, 
wc have Us = {2eB/h'^c) J^'dp^{exp[{pl/2me - fi)/kT] + l}'^ ~ {2eB/h'^c){2meHy/'^[l - 
(7rA;T)V(24//2)], which gives fj,c^Ep[l + {nkTf / {UEj,)]. Note that this correction has the 
same form as that for a non-magnetic electron gas, but with opposite sign. 



18 



where the cohesive energy Qg = \Esfi\ — |-E'(H)| is given by equation (2.16). Since PVsfi 
is typically much smaller than the cohesive energy Qs, and for low vapor density we can 
neglect the non-ideal gas correction A//(H), equation (4.4) reduces to an explicit expression 
for n(H): 

^^^^ " J [-wj '^^^ [ — T — J • ^^-^^ 

The densities of the other species in the saturated vapor can also be obtained. Equation 
(4.5) together with the Saha equation for e+p=H-|-(5i (Eq. [3.14]) yields 



rip ~ 



(27r)3/2 



tanh 



2MT 



-1/2 



exp 



-Qi - Qs + ^fj-s 

2T 



(4.6) 



The partition function of Hjv molecules can be approximated by 

3/2 



Z(Hjv) = V, 



(NMTy (M'^\ . , / NEn\ 



(cf. §3.3, and recall that En is the energy per atom in a Htv molecule). The equilibrium 
condition A^//s = f-i-N for the process Hg^oo + H = Hg^oo+Ar (with energy release Sn = 
NEn — NEg, where Sj\[ is the surface energy discussed in §2.5.) would then give n(Hjv) = 
exp(//Ar/r)Z(HAr)/Vg. For large molecular chains {N » iVj or droplets {N ^ Nc, cf. §2.5), 
however, a correction A/Xjv to the "internal" chemical potential due to finite temperature 
need to be included {Afi^ is to be distinguished from the non-ideal gas correction discussed 
in §3.4). For large 3d droplets, this correction is identical to that given in Eq. (4.3), i.e., 
AfiN = NA/ig- For linear chains with N » Ns, we expect that Ajjln /N is still close to AjjLs 
since electrons behave like a Fermi gas in the molecular chain as they do in 3d condensate, 
although it is not clearly which is greater.^ We therefore write A^^r = CNA^g, with ^ ~ 
for N ^ Ng and ~ 1 for >> A^^- With this correction, the saturated density of Hjv is 
given by 



-SN + N{l-OAfig 



For concreteness, we shall set ^ = for A'' = 1 — 5 and ^ = 1 for iV > 5 in our calculations 
below. This is adequate for our purpose since A'^(l — (^)A^s is typically much less than Sn 
for small molecules. More detailed (and unknown) prescription for ^ as a function of N (as 
long as ( satisfies the asymptotic behavior discussed above) would have negligible effect on 
our results presented below. 



4.2 Critical Temperature for Phase Separation 

The critical temperature Tcrit, below which phase separation between condensed metal 
and gaseous vapor occurs, is determined by the condition Ug = Ug = Up + n(H) + 2n(H2) + 
3n(H3) + • • •. In equation (2.17) for the surface energy, the factor / is of order unity, and 

^ Since the mean electron density in 3d condensate is higher than that in Id chain by 
about 0.2% (cf. Appendix A), we may expect that Ajjlm/N is larger than Ajig by ~ 0.4% 
according to Eq. (4.3). 
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approaches a constant for N » Ng. Substituting Eq. (2.17) with / ~ constant into Eq. (4.7) 
with the approximation 2;(Hjv) ~ z{B.), we obtain 



n.=^n,+ (^^J (^-^J-(H)exp 



/Qoo^ ^^5/2 ^3^pf iV[aQoo - (1 - QA/x.] j 

(4.8 



Clearly, for the sum to converge, we require aQ^o > (1 — C)A^s. With aQ^ = AEg — 
0.01|£'oo| ^ 2 5°2^7 eV, this reduces to T5 ^ 1.5 B^^^ /{I - C)^/^ which is easily satisfied for 
r < Tcrit (see below) and ~ 1. Thus we neglect (1 — C)A/Us in comparison to ckQoo- The sum 
in Eq. (4.8) therefore reduces to a finite number of order N^l"^ [exp {aQ^o/T) — l]""*^ ~ _/v5/2 
(recall that /Qoo >> ctQoo)-, where N ^ 2.5T/(a(5oo) specifies the term that contributes 
most to the sum. The critical temperature for phase transition is then given by 

r„,..Z^.„,l/<5^, A.„,-.(^)"^(^),(H,iV=/., (4..) 

where we have neglected Up, and in estimating In A we have used = Ps/nip ~ 50 5^2^ 
(a.u.) (cf. Eq. [3.1]) and ~ 2 3^^ ~ T„it (which gives TV ~ 20). Thus for / ~ 1, we find 
Tcpit ~ 10^ 5 X 10^ and 10^ K for B12 = 10, 100 and 500 respectively (see also Fig. 5). The 
corresponding critical pressure is of order Pcrit — {jig/ {N))kTcnt ^ O.lnsQoo/i^), where 
(N) is the typical size of molecules in the vapor. 

4.3 Saturated Vapor of Condensed Metallic Hydrogen 

Figure 5(a)-(c) depicts the partial saturation vapor densities of different species in equi- 
librium with metallic hydrogen as a function of temperature for B12 = 10, 100 and 500. These 
densities are calculated using Eqs. (4.5)-(4.7), where the surface energy Sn = {Na + f)Qoo 
is obtained from and Qoo via Eq. (2.18), and we adopt the same ansatz for Qn as de- 
scribed in §3.5, i.e., we use the numerical results for Qi, Q2, Qs, Qa and Qoo (the zero-point 
energy corrections to Q3 and Q4 are incorporated via Qn /Qoo = Qn°^ /Q''^^)i and we set 
Qn = Qa for A'^ > 4. This approximation for Qiv corresponds to / = 0.23, 0.6, 1.1, 1.7 for 
A?" > 4 at B12 = 1, 10, 100, 500 respectively. In calculating the total baryon density Ug in 
the vapor, we include linear molecular chains Hjv with A^ up to N^ax = 200 [cf. Eq. (2.19)] 
and neglect "droplets" in the vapor. Choosing a different N^ax >> 1 would only change the 
total vapor density Ug slightly for T <, Tcrit- The solid vertical line in each panel of Fig. 5 
indicates the temperature at which Ug equals n^. This defines the critical temperature Tcrit 
for phase transition, and it is well approximated by 0.1 Qoo/k. Thus Tcrit — 9 x 10^, 5 x 10^ 
and 1.4 x 10^ K for B12 = 10, 100 and 500 respectively. Note that for illustrative purpose, 
we have considered in Fig. 5(a) relatively low field strength, B12 = 10, where the condensed 
phase is more uncertain. This uncertainty for the low-field case is reflected by the fact that 
n(H2) > n(H) in the saturated vapor, so that larger molecules have even greater abundances. 
Thus there may not be a distinct condensed phase at all, and even if phase separation exists, 
the condensation temperature is below 10^ K. In such relatively low-B regime, the outer 
layer of the neutron star is characterized by gradual transformation from nondcgcnerate gas 
to degenerate plasma as the pressure (or column density) increases (sec §3). For B12 >> 10 
[cf. Fig. 5(b)-(c)], on the other hand, the vapor density becomes much less than the conden- 
sation density as the temperature decreases and/or the magnetic field increases, thus phase 
separation is inevitable. 



20 



The most abundant species in saturated vapor at T ~ Tcrit/'^ (for example) depend on 
the field strength, with smaller fields favoring poly-atomic molecules since / in Eq. (2.17) 
is smaller for smaller B. The ratio n(H2)/n(H) depends exponentially on [Si — S2)/T = 
(Q2 - Qs)/T ~ (Q2 - Qoo)/T. Prom Figure 1, we see that Qoo - Q2 > for for S12 ^ 10, 
and increases with increasing B. Therefore, when B12 » 10, there are few H2 molecules 
compared with H atoms in the saturated vapor. The ionization ratio np/n(H) in the vapor 
depends exponentially on [Qg — Qi)/{2T). Prom Figure 1, we see that Q^o — Qi > when 
B\2 ^ 200, in which case the vapor mostly consists of ionized hydrogen (cf. Fig. 5(c)). For 
10 ^ B12 ^ 200, n(H) is greather than Hp and n(H2), although the abundances of larger 
molecules are also appreciable. 

The column density j/sat above the surface of the condensed phase is related to the 
saturation vapor pressure Pgat by ysat = Psat/g and is plotted against T in Figure 6 for 
B12 = 10, 100, 500 (all for §14^ = 1). Note that T is the temperature at the gas/metal phase 
boundary. The density of the vapor is typically p ~ mpgy/{kT) ~ K^gnT^^y (g cm~^). 
Using equation (3.6) for the magnetic free-free opacity, we obtain the optical depth of the 
vapor above the metallic hydrogen, Tff ~ 4QQI3gii{T^^^^ B^2y^ ■ The threshold vapor column 
density yth, below which the vapor is optically thin, is then given by 

yth ~ 0.05r'/V4'^'^i2T5'/' (gcm-3). (4.10) 

This threshold value yth is also plotted against T in Fig. 6. The intersect of the ysat{T) and 
yth{T) curves defines a "threshold photospheric temperature" rph,th for the optical depth 
above the the surface of the condensate (if any): For Tph > Tph,th, the photosphere is purely 
in the gaseous phase, as described in §3. For Tph < Tph.th the vapor above the condensed 
metal is optically thin, and there is a sharp change in density at the gas/metal interface; The 
"photosphere" is located inside the condensed phase. For B12 = 500, the saturated vapor 
is dominated by ionized hydrogen, hence the vahie for Tph,th should be fairly accurate. For 
B12 = 10 and 100, however, the vapor is dominated by atoms or molecules, free-free opacity 
is only an estimate and thus the values of Tph,th given in Fig. 6 are more approximate. In 
any case, we see that rph,th is only slightly less than the condensation temperature Tcrit • 

In the case of Tph < Tph,th, the pressure at the gas/metal interface Pgat is less than Pph. 
Some distance into the condensed phase, with pressure increasing from Pgat to Pph (which 
is still much less than P^rit for the condensed phase), the liquid is still optically thin and T 
stays close to Tph. Further still into the metal, as P increases toward Pcrit from the "liquid 
photosphere" , the temperature increases appreciably while the density increases only slowly. 
Whether this region is convective or radiative depends on the opacity law in the condensed 
metal, which we have not yet investigated. 

5. DISCUSSIONS 

Recent works on the equation of state of neutron star surface in strong magnetic field 
have focused on the outer crust consisting of iron-like elements (e.g., Fushiki et al. 1989;Abra- 
hams Sz Shapiro 1991). Heat transport through this crustal region can be affected by the 
strong magnetic field (e.g., Hernquist 1985). However, for a neutron star covered by hydrogen 
atmospheric layer, the outer crust of the solid lattice is not the directly observable region 
of the star. Recent theoretical modeling of neutron star atmosphere aims at interpreting 
ROSAT's soft X-ray spectra from several radio pulsars, and most studies have focused on 
ionized atmospheres (Pavlov et al. 1995; Zavlin et al. 1995). 
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Our present study indicates that for sufficiently low temperature (T ^ 10^ K) and/or 
high magnetic field {B ^ 10^^ G), hydrogen atoms or molecules can have large abundance in 
the photosphere (see Fig. 2). This relatively large neutral abundance (even at relatively high 
temperature) comes about mainly for two reasons: (i) The binding energies of the bound 
states are greatly increased by the magnetic field; (ii) The photospheric density is much larger 
because the surface gravity is strong and also because the opacity of one of the two photon 
modes is significantly reduced by the magnetic field. Therefore one could in principle expect 
some atomic or molecular line features in the soft X-ray or UV spectra. In particular, the 
Lyman ionization edge (160 eV at 10^^ G and 310 eV at 10"^^ G) is expected to be prominent 
by the following consideration: The free-free and bound-free cross-sections for a photon in 
the extraordinary mode (with the photon electric field perpendicular to the magnetic field) 
are approximately given by (Gnedin, Pavlov &; Tsygan 1974; Ventura et al. 1992; Potekhin 
& Pavlov 1993): 

/Q \3/2 

where a = 1/137 is the fine structure constant, Oq is the Bohr radius, p is the density in 
gcm~^, Lo is the photon energy in the atomic units, Qi is the ionization potential, and b is 
the dimensionless field strength defined in Eq. (1.1). Near the absorption edge, the ratio of 
the free-free and bound-free opacities is 

^ ^ 10-4 P 



Thus even for relatively small neutral H abundance ^(H), the discontinuity of the total 
opacity at the Lyman edge is pronounced.^ Clearly, the position and the strength of the 
absorption edge can provide a useful diagnostic of the neutron star surfaces. 

For even stronger magnetic field {B larger than a few times lO"*^'^ G), we have shown that 
the degenerate matter can extend all the way to the outer edge of the star with little gaseous 
atmosphere above it, if the effective surface temperature is less than a threshold value 
Tph,th (cf. Fig. 6). In this case, the thermal radiation arises from just below the metallic 
hydrogen surface. Using equation (3.1) as an estimate for the condensation density (near 

zero presure), we find that the electron plasma frequency is given by hQ.e ~ 0.66 i?^2^ 
keV, much larger than the typical photon energy kT c± 10 Ts eV. Thus photons cannot 
be easily excited thermally inside the metallic sTirface. As a result, the surface emission is 
reduced from the black-body emission with the same temperature, and the spectrum of the 
surface radiation is likely to deviate quite strongly from a black-body spectrum (e.g., Itoh 
1975; Brinkmann 1980). We have not investigated the opacity structure of the hydrogen 
condensate and cannot say whether the emission mimics highly diluted blackbody radiation 
with color temperature much larger than the effective surface temperature. 

Finally, we note that the possibility that the outmost layer of a neutron star consists 
of degenerate iron has been considered previously, based on the notion that in a strong 
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Since the extraordinary mode has smaller opacity, most of the X-ray flux will come out in 
this mode. For photons in the ordinary mode, we have CTfj y c± 1.65 x 10^ pT^^^^ a'^ aluj~'^ (the 
zero-field value), and crbfy ^ 10^7rQ;ao(Qi/a;)^/^(ln6)-^, which give Kffj|/Kbfll ~ 10~'^ pT^^^^ 
for Bi2 ~ 1. Thus the absorption edge for the ordinary mode is less pronounced. 



22 



magnetic field the cohesive energy of Fe solid can be as large as tens of keV (Ruderman 
1974). While it is certainly correct that a strong magnetic field can significantly enhance 
the cohesive energy, it has been shown (Neuhauser et al. 1987) that the earlier calculation 
(Flowers et al. 1977) greatly overestimated the cohesive energy of Fe solid (cf. §2.3.2). Thus 
the metallic iron surface scenario was replaced by a non-degenerate atmosphere model. Our 
results in this paper partially resurrects the "degenerate surface picture": hydrogen (not 
iron) can condensate at sufficiently high B field and low temperature. 

This work has been supported in part by NASA Grant NAGW-2394 and NAG 5-2756 

to Caltech, and by NSF Grant AST 91-19475, NASA Grant NAGW-666 and NAG 5-3097 
to Cornell University. D.L. also acknowledges the support of a Richard C. Tolman Research 
Fellowship in theoretical astrophysics at Caltech. 



APPENDIX A. ESTIMATE THE RELATIVE BINDING ENERGY 
BETWEEN ID CHAIN AND 3D CONDENSATE 

In the absence of self-consistent electronic structure calculation for the 3d condensed 
metal, we estimate the fractional relative binding energy, A.Es/\E^\ = (Eg — £'(^)/|£'(^|, 
using the following three methods: 

(i) Method 1: If we approximate the Id chain by an uniform cylinder (radius R, proton 
spacing a), then the energy per atom (subcylinder) in the chain can be written as 

E =^ IL^-^-^^ 

°° 362i?4a2 a[ R y Aj 

where 7 = 0.57721566 (see Paper I). Minimizing E^o with respect to R and a gives 

R = 1.694 a/R = 1.885, E^ = -0.3913 6^/^ (A2) 

Comparing with equation (2.8) for the energy of 3d metal in the Wigner-Seitz approximation, 
we find AEJlE^ol ^ -1%. 

(ii) Method 2: Imagine that the 3d condensate is formed by placing a pile of parallel 
chains in contact with each other. The chain is assumed to have the property given by equa- 
tion (A2). For a body-centered tetragonal lattice, the relative position of nearest neighboring 
protons is xq = 2R perpendicular to the field and Zq = a/2 parallel to the field. We calculate 
the electrostatic interaction energy £"12 between a pair of nearest neighboring subcylinders 
using a Monte-Carlo integration method (Interaction between subcylinders of larger separa- 
tions can be calculated using a quadrupole-quadrupole formula, but this interaction is small 
compared to the nearest neighbor interaction). This gives £'12 = —3.7 x 10~'*6^/^. Since each 
proton has eight nearest neighbors, we have AEg = 4£i2, and thus AEs/\Eoo \ — —0.4%. 

If we choose xq < 2R (which would involve overlap of electron clouds), the chain-chain 
interaction is stronger. For xq = V2R (the smallest value possible), we find AEs/\Eoo\ — 
—8%. We consider this as an upper limit to the true relative binding energy. 

(ii) Method 3: Similar to (ii), we form 3d solid by placing chains together into a body- 
centered tetragonal lattice, with xq = 2R and zq = a/2. But instead of equation (A2), we 
determine R and a by minimizing the total energy of the solid. Eg = E^o + AEg, where E^o 



(Al) 
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is given by equation (Al), and the numerical result for AEg = AE12 can be fitted by the 
following analytic formula (with accuracy <, 1% for 1.6 < a/R< 2.2): 



AE"^ = -4.0 X 10"'* 4 
R 



1 + 23 



Minimizing Es{R,a), we obtain: 

R = 1.650 6-2/^ a/R = 2.036, E^ = -0.3932 6^/5. 
Comparing with equation (A2) for a single chain, we find AEs/\Eoo\ — —0.5%. 



{A3) 



APPENDIX B. PYCNONUCLEAR REACTIONS INDUCED BY 
STRONG MAGNETIC FIELD 



In this appendix, we consider the interesting possibility of "cold fusion" of hydrogen 
bound states induced by strong magnetic field on the neutron star surfaces. The magnetized 
H molecules and condensed metal arc highly compressed, i.e., the internal electron density 
within these bound states is large. This results in strong screening to the ion-ion Coulomb 
potential, and, as we will show, dramatic increase in the fusion rate between the nuclei within 
the bound states. 

Nuclear reaction in bound molecular systems has been studied in the context of muon- 
catalysed fusion, where a massive muon replaces the electron in a diatomic molecule of hy- 
drogen isotope, enhancing the binding and producing large cold fusion rates (e.g., Zel'dolvich 
& Gcrshtein 1961). In astrophysics, nuclear reactions induced by high matter density ('py- 
cnonuclear reactions') have also been studied (e.g., Salpeter & Van Horn 1969; Schramm 
&; Koonin 1990). The magnetic field induced fusion is similar to the density-induced pyc- 
nonuclear reaction except that the large density increase comes from the formation of bound 
states in strong magnetic field. 

We consider the enhancement of the fusion rate by magnetic field in both molecules and 
condensed metals. 



B.l Internal Fusion of H2 

Consider a diatomic molecule composed of two isotopic nuclei of hydrogen. The fusion 
rate A is proportional to the probability density |^'(r„)p that the two nuclei are at contact, 
i.e., A = A|\E'(r„)p, where A is the nuclear reaction constant which is directly related to the 
astrophysical S factor, and r„ is the nuclear radius. The relative wavefunction ^ of the two 
nuclei is governed by the interatomic potential V{r), which has been calculated in Paper I 
and II in the strong magnetic field regime. We can write V(r) in the form: 

V{r) = -- — + E{R2)-AVir), (Bl) 

r ro 

where tq is the equilibrium separation, E(H2) is the equilibrium energy of the molecule, and 
the function AV{r) = at r = ro and increases to [E(R2) — 1/ro — E(Ke)] as r decreases to 
zero {E{He) is the ground-state energy of He atom). For the molecule in the ground state, 
the relative energy Er of the two nuclei is simply E{}i2) plus the zero-point energy huj\\/2 
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(we consider vibration along the field axis only). In the WKB approximation, |\E'(r„)p is 
proportional to the WKB penetration factor P, i.e., 

|*(rn)P oc P = exp(-T^), W = 2^/2^1 j \v{r) - Er^/^dr, {B2) 

where Tc is the classical inner turning point, n is the reduced mass of the two nuclei in units 
of the electron mass. 

Neglecting Ay(r), equation (B2) can be easily integrated to give W ~ i^i^nr^^l"^ ^ with 
Tc = (l/^'o + ^i^ll/^)"^- A fitting formula for the numerical values of rg given in Paper I and 
Paper II is rp 12.7(ln6)~^'^. With ^ tq we have W ~ 15.8-y/7I(ln fe)"-*^'-^. The numerical 
results calculated using the exact potential Vir) and zero-point energy (see Eq. [2.4]) are 
given in Table 2, where the values of A = W/y^ for different field strength are listed (A does 
not depend sensitively on //; e.g., for Bx2 = 1, A ~ 1.74 for pp molecule and A ~ 1.75 for dd 
molecule.) These numerical results can be fitted by 

W = 15.3V^(ln6)-^■^ (B3) 

which is reasonably close to the approximate analytic expression. The factor A and the fusion 
rates for zero-field molecules as calculated by Koonin & Nauenberg (1989) are also listed in 
Table 2 for comparison. The reactions considered are: p-|-p^^H-|-e"'" -\- fg, p-l-d^'^He-l-7, 
p-|-t— >'^He-|-7, d-|-d^^He-|-n©^H-|-p, d-|-t^^He-|-n. Rescale the zero-field fusion rates using 
the appropriate A for strong field, we immediately infer the fusion rates for molecules in strong 
magnetic field. ^ These are listed in Table 2. Clearly, the strong magnetic field increases the 
fusion rate by many orders of magnitude. The energy generation rate per gram of molecules 
is ~ lO^^A ergg~^s~^. 

B.2 Internal Fusion of H^r Chains 

Larger Hjv molecules (before saturation; see §2.2) have similar electronic potential and 
excitation energies as H2. Therefore similar enhancement to the fusion rates also occurs in 
Hjv molecules, and this enhancement becomes more pronounced as N (and the field strength) 
increases. For a long chain molecule Hjv with 1 << N « Ng ~ [6/(ln ft)^]-*^/^, the spacing 
To along a field line between adjacent protons decreases with increasing N approximately as^ 
ro ~ 1/{N'^ In 6) (Paper I), and thus W ~ r^"^ ~ 1/N. A deviation 5r from the equilibrium 
spacing would give an excess potential 5V ~ (ln5)((5r)^/rQ. Thus the fractional zero-point 
vibration amplitude Ar/rQ is of order {rrif./mpY^^N^^'^ , i.e., the aligned vibrations become 
more pronounced as A'^ increases. Both the decrease in ro and the increase in Ar/ro tend to 
enhance the nuclear reaction rate. 

B.3 Fusion in Metallic Hydrogen 

As mentioned in §2.3, the zero-pressure metallic hydrogen has nonuniform electron dis- 
tribution. Although this screening effect will certainly increase the fusion rate, the essential 

^ In the strong magnetic field, the wavefunction \I' for the two nuclei is not spherical, 
unlike the filed-free case. Therefore the proportional factor in A oc ex.p{—Xy/]I) is different 
for B = and for high field cases, and our values for the fusion rate at high B are not exact. 
However, since the dominant factor determining the fusion rate is still the penetration factor, 
the error should not exceed an order of magnitude. 

^ Note that this is the asymptotic limiting result. For small N, the scaling with b is closer 
to ro ~ {\nb)~^. 



25 



TABLE 2 



Bound state fusion rates for isotopic hydrogen molecules 





Reactions " 


ir ir 


pd 


pt 


dd 


dt 




11 /nip 


1/2 


2/3 


3/4 


1 


6/5 




logio ^ (cm^s~^) 


-39.1 


-21.3 


-20.3 


-15.8 


-13.9 










logio A s"^ 









4.13 


-64.4 


-55.0 


-57.8 


-63.5 


-68.9 


1 


1.74 


-32.9 


-18.7 


-19.3 


-19.0 


-20.2 


5 


1.35 


-27.8 


-12.7 


-13.0 


-11.8 


-12.2 


10 


1.22 


-26.1 


-10.8 


-10.9 


-9.3 


-9.6 


100 


.895 


-21.8 


-5.8 


-5.7 


-3.3 


-3.0 


500 


.740 


-19.8 


-3.5 


-3.2 


-0.4 


0.2 



" The reactions are given in the text. 

A is defined via P = exp(— Ay^); and the values are given for the pd system; see 
text. 



feature can be obtained in the uniform electron gas model. As an estimate, we can simply 
use the fitting formula of Salpeter and van Horn (1969; sec also Schramm & Koonin 1990) 
for the pycnonuclear reaction rate, except that the zero-pressure density is given by equa- 
tion (3.1). For pd reaction, the rate per deuteron is Apj ~ 2 x 10^^pg^^^exp(— 21.8/9q ^''®) 
s~^. Thus a deuteron embedded in a cold hydrogen plasma has lifetime A~J = 10^-^ 
years at B12 = 2 and one year at B12 = 13. For pp fusion, the reaction rate is App c±. 
800 /9g^^^exp(— 18.89/9g "^^^) yr~^. At B12 = 60, 100, 500, corresponding to zero-pressure 
densities p<,,o - 7.6 x 10"^, 1.4 x 10^ 10^ g cm'^, we find A'p^ ~ 10^^-^, 10^, lO^-^ yrs, re- 
spectively. Thus for non-accreting neutron star, there is an upper limit to the field strength 
above which the surface hydrogen layer is absent. Comparing with Table 2, we see that 
at a given field strength, the fusion rate is larger in the metallic state than in molecule. 
The reason is that the mean ion separation in the condensed state decreases with increasing 
field strength as a power-law function of B, whereas the mean interatomic separation in a 
molecule decreases only logarithmically. Also note that as the matter density or temperature 
increase, the magnetic field effects become less important because many Landau levels are 
occupied by the electrons. 
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FIGURE CAPTIONS 



FIG. 1. — Energy releases from several atomic and molecular processes as a function of the 
magnetic field strength. The solid line shows the ionization energy Qi of H atom, the dotted 
line shows the dissociation energy Q2 of H2, and the dashed line shows the cohesive energy 
Qoo of linear chain Hoo- The zero-point energy corrections have been included in Q2 and 

FIG. 2. — Temperature-density diagram of the atmosphere of a neutron star with sur- 
face field strength (a) B12 = 1 and (b) B12 = 10, both with 514 = 1. The solid lines 
correspond to the photosphere ("r ~ 2/3") and the extraordinary photon emission region 
("r_L ~ 2/3"); the dotted curves ("e+p=H") correspond to constant neutral H fractions, 
with n(H)/ng = 0.9 for the lower curves, n(H)/ng = 0.5 (half/half ionization) for the middle 
curves, and n{H)/ng = 0.1 for the upper curves; the dashed curves ("H-|-H=H2") correspond 
to the half/half dissociation of H2; the long-dashed lines ("T = Tp") show the Fermi temper- 
ature; the dot-long-dashed lines show selected temperature-density profiles of atmospheres 
with different photospheric temperatures (Tph,5 = 5, 2 for B12 = 1 and Tph,5 = 9, 3 for 
B12 = 10); the dot-dashed vertical lines ("p = PLandau") correspond to the density above 
which the electrons occupy the excited Landau level; the dark vertical lines correspond to 
the condensation density o! The filled circles correspond to the photospheric "ionization 
temperature" Tph,ir (see text). 

FIG. 3. — Neutral fraction n(R)/ng as a function of temperature at a fixed density p = 
HgTUp = 0.01 g/cm^ (heavy lines) and 0.1 g/cm^ (light lines). The curves are labeled by the 
field strength: B12 = 0.01 (solid lines), 0.1 (short-dashed lines), 1 (long-dashed lines), 10 
(dot-dashed lines). The zero-field results are also shown for comparison (dotted lines). 

FIG. 4. — Abundance of various H species as a function of the photosphere temperature 
at a point with P = Pph/2 (just a little above the photosphere) for three field strengths 
(a) B12 = 1, (b) 10, and (c) 100. The fractional number density are defined by A(Hjv) = 
n(Hiv)/ntot5 = np/ntot, with ntot = 2np -I- n(H) -|- n(H2) -I- • • •. The numbers labeling 
each curve specify N, with, e.g., "3-5" indicating ^(Ha) -|- X{B.4) + ^(Hg). 

FIG. 5. — The saturation vapor densities of various species (in the atomic units, a~^) 
of condensed metallic hydrogen as a function of temperature for different magnetic field 
strengths: (a) B12 = 10; (b) B12 = 100; (c) B12 = 500. The dotted curves give Hp, the 
short-dashed curves give n(H), the long-dashed curves give n(H2), the dot-dashed curves 
give [3n(H3) -|- 4n(H4) -I- • • • -I- 8n(H8)], and the solid curves give the total baryon number 
density in the vapor Ug = Up + n(H) -|- 2n(H2) + • • •. The horizontal solid lines denote the 

condensation density ~ 50 5^2^ (a.u.), while the vertical solid lines correspond to the 
critical condensation temperature at which Ug = Ug- 

FIG. 6. — Column density of saturated vapor. The heavy (steeper) lines show the column 
density y of the nondegenerate gas above the condensed metallic phase as a function of the 

temperature at the phase boundary, the lighter lines show the threshold column densities 
Uth (Eq. [4.10]) below which the vapor is optically thin to free-free absorption. The solid 
lines are for B12 = 10, the short-dashed lines for B12 = 100, and the long-dashed lines for 
B12 = 500, all with 514 = 1. The filled circles and vertical dotted lines correspond to the 
threshold photospheric temperature Tph^th below which the saturated vapor is optically thin. 
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